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ABSTRACT 

We explore here the possibly important influence of cooled interstellar gas on the fundamental 
plane of elliptical galaxies. Interstellar cooling is described by a parameterized sink term in the 
equation of continuity. Parameters that give the best fits to the X-ray observations of elliptical 
galaxy NGC 4472 are used as a template for other galaxies. These spatially-dependent cooling 
parameters are then extended homologously to elliptical galaxies of lower mass for which X-ray 
observations are currently unavailable or impossible to observe because of stronger relative X-ray 
emission from binary stars. Most of the interstellar gas cools within an effective radius where it 
can contribute an additional 10 - 30 percent to the mass of the old stellar population. The absence 
of observed cooled gas and simple theoretical arguments support the hypothesis that the cooled 
gas is forming into low mass stars, with implications that depend on the IMF of this young stellar 
population. If the IMF includes only stars of very low mass, <C Mq as is commonly believed, 
the cooled mass is optically dark. For this IMF the mass to light ratios determined from stellar 
velocities systematically overestimate that of the old stellar population. Furthermore, the total 
mass and spatial distribution of the optically dark stellar population does not scale homologously 
with galactic luminosity or radius so the total stellar mass to light ratio is expected to vary 
with galactic radius. These variations in stellar population and non- homology can introduce 
pronounced deviations from the fundamental plane. We investigate cooled gas perturbations to 
the mass to light ratio for several idealized homologous elliptical galaxies and show that they 
may be incompatible with the observed thinness of the fundamental plane. However, if the IMF 
of the stellar population produced from cooled interstellar gas is optically luminous, its influence 
on observed stellar mass to light ratios and the fundamental plane would be lessened. If the IMF 
of a young, continuously forming and optically luminous stellar population were sufficiently well- 
behaved, it could account for about ten percent of the optical light from large elliptical galaxies 
within r e , possibly having important implications for understanding the fundamental plane. 

Subject headings: galaxies: elliptical and lenticular - galaxies: structure - galaxies: fundamental plane - 
galaxies: cooling flows - x-rays: galaxies 



1. INTRODUCTION 

The mass of hot interstellar gas that cools over 
a Hubble time in a massive elliptical galaxy is typ- 
ically 5-7 percent of the total baryonic mass of 
these galaxies. Several observations indicate that 
this gas cools in a significant volume of the in- 



ner galaxy: (1) the total cooled mass significantly 
exceeds the mass of central massive black holes 
(Magorrian et al. 1998), (2) although much of the 
interstellar gas within r e is ejected from old galac- 
tic stars which typically have a systematic galactic 
rotation, the inner X-ray images of rotating ellipti- 



cal galaxies are nearly circular, implying that the 
gas cools before it flows very far toward the ro- 
tation axis (Hanlan & Bregman 2000; Brighcnti 
& Mathews 2000b), (3) direct evidence of diffuse 
cooled gas at T ~ 10 4 K is invariably observed 
in large elliptical galaxies in r ^ 0.3r e (e.g. Mac- 
chetto et al. 1996), and (4) evolutionary models 
for the hot interstellar gas are in better agreement 
with X-ray observations when cooling mass loss is 
included; most of this gas cools within an effective 
radius (Brighenti & Mathews 2000a). 

In view of the very small masses of cold HI or H 2 
gas that have been observed in elliptical galaxies, 
far less than the mass of hot gas that has cooled, it 
is usually assumed that the cooled interstellar gas 
is continuously forming into a population of very 
low mass (brown dwarf) stars (Fabian, Nulsen & 
Canizares 1982). Such stars add mass but no light. 
We show here that if the cooled gas is optically 
dark, the fundamental plane for elliptical galaxies 
is likely to be perturbed beyond observed limits. 
Alternatively, if the mass to light ratio of the form- 
ing stars is not infinite but is comparable to that 
of the old stars, then some of the cosmic scatter or 
tilt of the fundamental plane from the virial plane 
could be due to the presence of this young stellar 
population. Support for a continuously forming, 
optically luminous stellar population in elliptical 
galaxies is provided by the nearly solar mass of 
neutral cores in interstellar clouds at which gravi- 
tational collapse first occurs, ~ 2 Mq (Mathews & 
Brighcnti 1999a), and by the apparent young stel- 
lar ages of many elliptical galaxies based on the 
H/3 photometric index (Worthey 1994; Worthcy, 
Trager & Faber 1996; Trager 1997; Mathews & 
Brighenti 1999b). 

In the following discussion, however, we adopt 
the usual hypothesis that newly forming stars in 
elliptical galaxies are optically dark. In this case 
radial variations in the total stellar mass to light 
ratio must occur. Furthermore, the total mass of 
cooled gas varies with galactic luminosity, result- 
ing in rather profound excursions from the funda- 
mental plane, evidently in violation of the remark- 
able thinness and regularity of this plane. 

2. THE FUNDAMENTAL PLANE 

If elliptical galaxies were perfectly homolo- 
gous stellar systems with identical stellar popu- 



lations, then the virial theorem a 2 — k v (M/r e ) = 
n v {Ly /r e )(M / Ly) would imply simple correla- 
tions among observable parameters. Here a is the 
"central" velocity dispersion, M is the stellar (or 
galactic) mass, r e is the effective (half-light) radius 
and k v is a constant of proportionality. The three 
observables — cr 2 , Ly, and r e - define a homologous 
plane (HP) in this (logarithmic) three-dimensional 
parameter space. Instead, elliptical galaxies sat- 
isfy a different, slightly steeper fundamental plane 
(FP) relation a 2 cx (Ly /r e )[r°.- 22 a® A9 }, implying 
that k v (M/L v ) oc r e 22 a° o A9 oc M°- 24 r~ 002 cx 
L v °- 32 r-°- o:i (Dressier et al. 1987; Djorgovski & 
Davis 1987). In view of the importance of ellipti- 
cal galaxies for D n — a D distance determinations, 
considerable effort has gone into understanding 
the physical origins of the "tilt" of the FP relative 
to the HP. 

Interpretations of this tilt are constrained by 
the regularity of the FP-HP deviation across the 
entire region observed and the remarkably small 
observational scatter (~ 12%) about the FP rela- 
tion (Rcnzini & Ciotti 1993). Generally, the FP- 
HP tilt has been interpreted as some combination 
of two limiting effects: either (1) stellar popula- 
tions or dark matter contributions vary systemat- 
ically across the FP so that M/Ly oc Ly ' 32 or 
(2) homology fails so that k v oc Ly . This re- 
lation depends mildly on photometric band (e.g. 
Scodeggio et al. 1998; Pahre & Djorgovski 1997). 

Several types of stellar population variation 
(age, metallicity or IMF) have been considered 
along the FP (Renzini & Ciotti 1993; Pahre, 
Djorgovski & Carvalho 1995; Pahre & Djorgov- 
ski 1997). These attempts to interpret the FP- 
HP tilt have not been entirely successful since the 
FP is nearly invariant to observed wavelength over 
the B — H photometric bands. The small system- 
atic wavelength trend of the FP slope is consistent 
with the color-magnitude relation for giant ellipti- 
cal galaxies (Scodeggio et al 1998; Pahre, de Car- 
valho & Djorgovski 1998). Moreover, if the stellar 
population changes across the FP, differential stel- 
lar evolution can in some cases introduce a change 
in the FP slope at large redshifts, but no such ef- 
fect has yet been observed. It is generally asserted 
that the full FP tilt cannot be understood solely 
by metallicity variations or changing stellar pop- 
ulations along the FP since this would require a 
high degree of fine tuning in the IMF. If M/L is 
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assumed to be constant with galactic radius, the 
mass to light ratios for elliptical galaxies found 
from dynamic models vary in a manner similar 
to the FP deviation, M/L oc iO.35io.05 ( van dcr 
Marel 1991). 

A wide variety of non-homologous structural 
and/or kinematic models have been proposed to 
explain all or part of the FP-HP tilt (e.g., Capelato 
et al. 1995; Djorgovski 1995; Hjorth & Madsen 
1995; Ciotti, Lanzoni & Rcnzini 1996). In these 
types of models the stellar MjLy docs not vary 
so the FP tilt should be strictly independent of 
wavelength apart from mutually non-homologous 
stellar and dark matter distributions. Graham & 
Colless (1997) fit the stellar mass and light dis- 
tribution with Sersic r 1 /" laws and note that n 
increases with galactic mass (by ~ 2 — 4) in a 
non-homologous fashion across the FP. If these 
galaxies are mapped onto the FP with an r 1 / 4 
law in the usual way, the resulting FP, r e oc 
(T i.33±o.io s -o.79±o.ii 5 tiltg even fui . thci . f rom the 

HP, r e oc o-qS^ 1 , where S e oc LJr\ is the surface 
brightness. However, a dynamical non-homology 
must accompany a structural non-homology. Solv- 
ing the Jeans equation for elliptical galaxies mod- 
eled with r 1 /" profiles (Ciotti 1991; Ciotti & 
Lanzoni 1997), Graham & Colless (1997) and 
Graham (1998) determine infinite-aperture Otot.n 
which when fit to r 1 /" galaxies result in r 6i „ oc 
^I 2 „ ±0 - 24 ^°- 74±0 - 09 , which is somewhat closer to 
the HP. Busarello et al. (1997) consider a different 
dynamic non-homology based on partial rotational 
support within r e . Using a modified a which incor- 
porates both rotation and random velocities, they 
find that more than half of the FP-HP tilt is ac- 
counted for, r e oc ^i-53±o.2 S -o.92±o.i Busarello et 
al. claim that the remaining discrepancy can be 
understood with spatial non-homology and stel- 
lar population differences (e.g. Prugniel & Simien 
1997). 

Another approach is the possibility that dark 
(halo) matter makes a variable contribution along 
the FP (Ciotti, Lanzoni, & Renzini 1996), result- 
ing in the FP-HP tilt. However, these authors 
stress again the tightness of observational scatter 
about the FP so that any systematic trend in the 
dark matter content must be rather finely tuned. 
This high degree of coordination may seem un- 
likely if elliptical galaxies are produced by a vari- 
ety of random mergers among stochastically differ- 



ent galaxies of lesser mass. The amount and radial 
distribution of dark halo matter in the centers of 
elliptical galaxies is complicated by the baryonic 
winds in first generation galaxies and the dissi- 
pational baryonic concentration expected as giant 
elliptical galaxies form by mergers; these counter- 
acting effects respectively contract and expand the 
dark matter distribution relative to the baryonic 
stellar population. 

There is at present no clear winner among the 
various explanations for the FP-HP tilt. Instead, 
a patchwork of small contributions - stellar pop- 
ulation variations, rotational support and non- 
homology - are required to bring FP into agree- 
ment with HP (Ciotti, Lanzoni, & Renzini 1996: 
Prugniel & Simien 1997; Pahre, de Carvalho & 
Djorgovski 1998). However, there is no compelling 
reason why these quite different and presumably 
independent contributions should combine in a co- 
ordinated fashion to produce the full FP-HP tilt 
observed across the entire FP. Moreover, a mul- 
tiplicity of coordinated explanations for the tilt, 
none entirely dominant, makes the tightness of the 
FP even more remarkable. 

3. COOLING DROPOUT MASS AND 
LOW MASS STARS 

In this paper we introduce a new complication 
to the FP-HP controversy: a variation in the dy- 
namical M/Lb with galactic mass due to the de- 
position of cooled interstellar gas in r < r e . The 
importance of cooled gas in the galactic mass bud- 
get follows from the universality of hot interstel- 
lar gas in giant elliptical galaxies and the direct 
relationship of the cooled mass to the loss of en- 
ergy implied by the observed X-ray emission L x . 
For example, interstellar gas in the large Virgo el- 
liptical galaxy NGC 4472 with temperature T w 
1.3 x 10 7 K emits L x (0.5 - 4.5kcV) = 4.5 x 10 41 
ergs s _1 in the Einstein band. The bolometric X- 
ray luminosity, L x ^oi ~ 1.6^^(0.5 — 4.5kcV), is 
the product of the mass cooling rate M and the 
specific enthalpy of the gas, therefore 

M=(^)i x , 6o; -2.5M yr- 1 . (1) 

Since L x was somewhat larger in the past, over 
a Hubble time - 3 - 4 x 10 10 M Q of hot gas will 
have cooled. Although this amounts to only about 
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~ 5 — 10 percent of the total stellar mass in NGC 
4472, M» t = 4.73 x 10 11 M , most of the cooled 
gas is likely to be deposited near the galactic cen- 
ter r < r e where the hot gas density and X-ray 
emissivity are greatest and where its contribution 
to the dynamical M/Lb determined from stellar 
velocity dispersions is most apparent. (In eval- 
uating M*t we assume M/Lb = 6, as explained 
below.) 

Although there is little doubt that interstellar 
cooling has occurred, two main questions must be 
addressed: (1) What is the final physical state of 
the cold gas? and (2) What is the radial dis- 
tribution of the cooled mass in the galaxy? For 
many years it has been assumed that low mass 
stars are the most likely final state for cooled gas 
(e.g. Fabian, Nulsen & Canizares 1982). In a re- 
cent paper (Mathews & Brighcnti 1999a) we de- 
scribe how the hot interstellar gas cools (initially) 
to T ~ 10 4 K where it can be heated and ion- 
ized by stellar UV photons and made visible in 
optical emission lines. This first stage of cooling 
occurs at many (~ 10 6 ) cooling sites distributed 
over the inner galaxy. However, when the radius 
of partially cooled gas at 10 4 K exceeds the local 
Stromgren radius at a cooling site, a cold neutral 
atomic or molecular core forms with T c ~ 15 K 
and n c ~ 10 5 cm~ 3 in pressure equilibrium with 
the ambient hot interstellar gas. As the mass of 
the neutral core M c increases further, it becomes 
Jeans unstable when M c ~ 1 — 2 M Q so only stars 
of this or lower mass can form from the available 
gas (Ferland, Fabian & Johnstone 1994; Math- 
ews & Brighenti 1999a). Cooling gas is processed 
within HII cooling site clouds with star-forming 
cold cores in time scales (~ 10 5 years) that are 
short compared to their hydrodynamic disruption 
times in the cooling flow atmosphere (Mathews & 
Brighenti 1999a). The high efficiency of low mass 
star formation nicely accounts for the absence of 
observable cold HI or molecular gas in large ellipti- 
cal galaxies, e.g. for NGC 4472 M HI + M H2 < 10 7 
Mq . We conclude that two stellar populations co- 
exist in giant elliptical galaxies: an old luminous 
stellar component with an IMF similar to Salpeter, 
and a second population of low mass stars contin- 
uously forming from the interstellar medium and 
having a bottom heavy IMF. Although the de- 
tailed character of this second IMF is unknown, for 
the purposes of the following discussion we adopt 



the usual assumption that it is restricted to very 
low masses, i.e. its mass to light ratio is essentially 
infinite. 

The second question - concerning the radial dis- 
tribution of cooled gas mass - can be answered in 
an approximate way by comparing hydrodynamic 
models with the radial properties of interstellar 
gas observed in X-rays. The production of low 
mass stars at many distributed sites can only occur 
if the interstellar gas is inhomogeneous. Radiative 
cooling is locally enhanced in regions of low en- 
tropy (low T, high density n). As gas in the galac- 
tic cooling flow slowly proceeds inward, regions of 
progressively higher entropy cool at smaller galac- 
tic radii; regions having low entropy cool further 
from the galactic centers. Distributed Ha + [Nil] 
line emission observed in the central regions of gi- 
ant elliptical galaxies illuminate the cooling sites 
where this process is occurring. Interstellar en- 
tropy variations can be created by a variety of 
complex events: stellar disturbances (e.g. stel- 
lar winds or Type la supernovae), interactions 
with merging dwarf galaxies or differential Type 
II supernova heating in the early universe during 
the earliest phases of star formation. Because of 
our poor understanding of these interstellar en- 
tropy fluctuations, it is not possible to determine 
from first principles where most of the hot gas 
cools and deposits its mass into low mass stars, 
although the observed optical line emission is a 
visible tracer of this process. The best currently 
available means of setting limits on the radial dis- 
tribution of cooling gas in elliptical galaxies is to 
compare the observed radial distribution of inter- 
stellar X-ray emission with the computed X-ray 
emission from gas-dynamic models for the cool- 
ing flow gas in which cooling dropout is assumed 
to occur according to various radial prescriptions. 
We have recently computed a series of models for 
the large elliptical galaxy NGC 4472 with a vari- 
ety of assumed mass dropout functions (Brighenti 
& Mathews 2000a). In the following section we 
describe the cooling dropout properties of models 
that best fit the X-ray observations, then we dis- 
cuss the implications of this dropout mass for the 
fundamental plane. 
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4. DESCRIPTION OF HYDRODYNAM- 
ICAL MODELS 

In view of the uncertainty concerning the cool- 
ing dropout, in our gas dynamical models of galac- 
tic cooling flows we parameterize the radial mass 
dropout of cooling gas by introducing an ad- 
justable sink term in the equation of continuity, 

M+^d- r ( r Pu) = ap* - q(r)-. (2) 

Here p(r, t) and u(r, t) are the gas density and ve- 
locity in the cooling flow. The sink term for cool- 
ing flow dropout is characterized by a dimension- 
less parameter q(r) and tdo — 5m p kT/2ppA is the 
time for gas at temperature T(r, t) to cool locally 
by radiative losses at constant pressure (see, e.g. 
Sarazin & Ashe 1989). The source term a(t)p*(r) 
represents the rate that interstellar gas is supplied 
by stellar mass loss from the old stellar population 
having density (r ). Since a(t) oc (M/Lb)^ 1 for 
a wide range of power law IMFs, the total stel- 
lar mass ejected J a(t)M*dt is nearly identical for 
elliptical galaxies with the same Lb but different 
IMFs (Brighcnti & Mathews 2000a). Following 
the approach in Brighcnti & Mathews (2000a), 
we consider galactic cooling flow models for NGC 
4472 using a parameterized dropout coefficient 

q(r) =q cxp(-r/r d0 ) m (3) 

where q a , r^o and m are adjustable parameters. 
For all sets of dropout parameters considered we 
found that the total dropout mass in r < r e /3 is 
a substantial fraction (0.10 - 0.35) of the mass of 
luminous stars in this region. The radial distribu- 
tion of the dropout mass, assumed to remain at 
the cooling site, does not follow the space density 
corresponding to the de Vaucouleurs light profile. 

If low mass stars formed from cooling dropout 
contribute mass but no appreciable light, the hy- 
pothesis that we want to test, it follows that the 
intrinsic "stellar" mass to light ratio of the lu- 
minous stellar population in NGC 4472 must be 
somewhat lower than values determined from the 
observed motion of luminous stars: M/Lb ~ 9.2 
(van der Marel 1991) or M/L B « 10.3 (Faber et 
al. 1997). The most satisfactory mass dropout 
coefficients q(r) generate hydrodynamic solutions 
for the interstellar gas at the current time (t n = 13 
Gyrs) that reproduce best the observed X-ray im- 
age of NGC 4472. Using this approach, we find 



that the dynamical mass to light ratio in successful 
hydrodynamic models requires a somewhat lower 
intrinsic mass to light ratio for the luminous stars 
M/Lb ~ 6. For the low mass dropout stellar pop- 
ulation we assume (M/LB)do — oo. 

Ideally, it would be desirable to determine if 
the same (or appropriately scaled) mass dropout 
coefficient q(r) used in gas dynamical models for 
elliptical galaxies of varying Lb, spread across 
the fundamental plane, can fit the observed X-ray 
surface brightness distributions T, X (R) in these 
same galaxies. Unfortunately, at the present 
pre-Chandra time sufficiently accurate H X (R) are 
known only for several bright elliptical galaxies 
in Virgo (e.g., Trinchieri, Fabbiano & Canizares 
1986), all having similar Lb- Approximate 
T, X (R) distributions are available only for about 
a dozen elliptical galaxies altogether (Mathews & 
Brighcnti 1998a). In addition, emission from low 
mass binary stars begins to dominate the X-ray lu- 
minosity in elliptical galaxies several times less lu- 
minous than NGC 4472, masking the X-ray emis- 
sion from the cooling flows. Finally, even though 
L x cx L 2 B approximately holds among brighter el- 
liptical galaxies, there is an enormous spread in 
values of L x for given Lb- We believe that most 
of this spread is due to environmental interac- 
tions: (i) tidal exchange of halo hot gas (and dark 
matter) in the group environments in which these 
elliptical galaxies formed (Mathews & Brighcnti 
1998b) or (ii) ram pressure stripping. 

With these difficulties in mind, we chose in- 
stead to explore the relative influence of optically 
dark interstellar dropout not directly on the FP 
but on a homologous galactic plane (HP) that is 
scaled from NGC 4472 for which good X-ray ob- 
servations are available. We consider a sequence 
of three massive elliptical galaxies, NGC 4472 and 
two similar model galaxies, 4472/4 and 4472/16, 
having optical luminosities Lb that are respec- 
tively 4 and 16 times less than NGC 4472. To 
explore only the structural non-homology caused 
by cooling dropout, we also assume that the old 
stellar population is invariant along our HP, i.e. its 
M/Lb is constant. For perfect homology we also 
require that the mass and scale radius of the dark 
halo distribution scale with the baryonic mass. Fi- 
nally, we solve for the evolution of hot gas in each 
galaxy along the HP using identical (or scaled) 
parameters (q , rdo and m) that give good agree- 
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ment with the X-ray observations for NGC 4472, 
T(r, t n ) and n(r, t n ); this procedure is described in 
more detail in Brighenti & Mathews (2000a). By 
this means we can determine the differential, non- 
homologous distribution of optically dark dropout 
mass and assess the deviation of these three ide- 
alized galaxies from the HP relation. We expect 
that these deviations from the HP arc similar to 
deviations of galaxies of similar Lb from the local 
FP. If our procedure were perfectly correct - i.e., 
elliptical galaxies are exactly homologous and con- 
tain identical old stellar populations - we would 
find that the non-homology introduced by cooling 
dropout would explain the "tilt" of the observed 
FP from the HP. Unfortunately, this goal has not 
yet been achieved. 

For perfect homology and simplicity we as- 
sume that 4472 and its two clones of lower mass 
are spherical with dc Vaucouleurs stellar profiles 
and scaled NFW dark halos (Navarro, Frcnk, & 
White 1996). The effective radius of NGC 4472 
is r e = 8.57 kpc (for distance d = 17 Mpc) and 
its total (old) stellar mass is M*_ t = 4.73 x 10 11 
M Q , based on M/L B = 6. The NFW halo in 
NGC 4472 has a virial mass of M h . vir = 4x 10 13 
Mq . These parameters give a good fit to the X-ray 
mass profile of NGC 4472 assuming hydrostatic 
equilibrium although the NFW profile is found 
to be somewhat too centrally peaked (Brighenti 
& Mathews 2000a). For the smaller comparison 
galaxies (4472/4 and 4472/16) we choose identical 
stellar parameters but scale the size of the stel- 
lar distribution downward by 4^ - 793 = 0.333 or 
^g-o.793 _ o.m where the exponent is found from 
the observed My-r e relation for massive elliptical 
galaxies, r e cx L°; 793 oc L° B 793 , (Faber et al. 1997). 
The virial mass of the NFW halos for the smaller 
galaxies are reduced by 4 or 16 and the scale radius 
for these halos varies as r s cx M^-^J , as required 
by our chosen cosmology for which the NFW con- 
centration varies as c oc M^ vir (Navarro, Frenk, 
& White 1996). 

As explained in more detail elsewhere (Brighenti 
& Mathews 1999a), our gas dynamical models be- 
gin with an overdensity in a simple Einstcin-dc 
Sitter cosmology (CI = 1, Clb = 0.05, and H a = 50 
km s _1 Mpc^ 1 ) toward which baryonic and dark 
matter are gravitationally attracted. Galactic 
stars are assumed to form at time t* s — 1 Gyr 
and, soon afterwards, Type II supernova energy 



based on a Salpetcr IMF (slope: x — 1.35, mass 
limits: mg = 0.08 and m u — 100 Mq) heats the 
ambient gas within the accretion shock with 80 
percent efficiency (see Brighenti & Mathews 1999a 
for details). At t* = 2 Gyrs the de Vaucouleurs 
stellar profile and potential is constructed, simu- 
lating the results of mergers. For this approximate 
exploration we ignore the possibility that ellipti- 
cal galaxies of lower luminosity may have formed 
slightly earlier. In some models we allow the sec- 
ondary inflow of cosmic gas to continue through- 
out the calculation. For comparison, in other 
models we ignore the secondary infall entirely, 
crudely simulating a tidal truncation and removal 
of halo material in a relatively dense group of 
newly formed galaxies, a process more likely in 
elliptical galaxies of moderate or low mass. 

5. RESULTS 

In Table 1 we list the main results of a vari- 
ety of hydrodynamic cooling flow models for 4472, 
4472/4 and 4472/16 with various assumed dropout 
parameters and assumed amounts of accumulated 
cosmic gas. Of particular interest are the mass 
to light ratios evaluated at r e /3 and at a fixed 
distance of 480 pc from the galactic centers. Re- 
call that M/Lb = 6 is assumed for the luminous 
stars in all three galaxies and that the observed 
value for NGC 4472 is - 9.2 (van der Marel 1991) 
or ~ 10.3 (Faber et al. 1997). Models in which 
the (stellar plus dropout) mass to light ratio for 
NGC 4472 falls in the range M/L B « 8.5 - 11 
can be regarded as reasonably successful. How- 
ever, the M/L B (r e /3) entries in Table 1 are en- 
hanced by about 15 percent due to a small amount 
of NFW dark halo material in this region. Since 
our fits to the X-ray surface brightness of NGC 
4472 are improved if the dark halo is less cen- 
trally peaked than NFW (also see Kravtsov et 
al. 1998; Burkert & Silk 1997; Navarro & Stein- 
metz 2000), we believe that its contribution to 
M/Lb in this region should be removed. There- 
fore, for comparison with Af/Ls(r e /3) values in 
Table 1 which contain this NFW excess, the range 
of acceptable models for NGC 4472 should be re- 
vised upward to M/L B (r e /3) « 9.8 - 12.6 (i.e. 
log[M/L B (r e /3)} = 0.99 - 1.10). Values of M/L B 
at 480 pc may be more relevant for comparison 
with M/Lb determinations that are often based 
on stellar velocity dispersions measured within a 
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fixed aperture of a few arcseconds about the galac- 
tic centers. 

In cooling flows without distributed dropout 
(q = 0) the flow proceeds entirely to the origin 
before cooling and the mass to light ratios at 480 
pc greatly exceed observed values for 4472 - such 
models can certainly be rejected. Evolutionary gas 
dynamic models with uniform q a — 1 at all radii 
and with dropout parameters (q , rdo, m) = (4, 2 
kpc, 1) both fit the observed distribution of inter- 
stellar gas temperature and density in NGC 4472 
reasonably well at the present time (t n = 13 Gyrs) 
as illustrated in Figure 1. While the agreement in 
Figure 1 is acceptable, it is by no means perfect. 
The apparent gas density (enhanced by dropout 
emission) exceeds observed densities in the central 
regions (r < 3 kpc) of NGC 4472, suggesting that 
some additional non-thermal support is present; 
NGC 4472 has a weak extended radio source in 
this region with an appreciable equipartition pres- 
sure (Ekers & Kolanyi 1978). At large galactic 
radii (r > 10 kpc) our models in Figure 1 are 
slightly lower than the azimuthally-averaged val- 
ues observed in NGC 4472; this density mismatch 
can be corrected with small changes in the galaxy 
formation time t* or the efficiency of SNII energy 
release (Brighenti & Mathews 1999a). However, 
our interest here is not to achieve absolute per- 
fection in fitting the observations of NGC 4472, 
but to make a differential comparison of 4472 with 
its less massive, homologous clones; 4472/4 and 
4472/16. The additional model with stronger but 
uniform dropout q = q a = 4, also included in Table 
1, results in apparent gas temperatures at t n = 13 
Gys that are considerably lower than those ob- 
served (Brighenti & Mathews 2000a); this is due 
to enhanced X-ray emission from cooling regions 
throughout the cooling flow. 

Several interesting details are apparent in Ta- 
ble 1. Since the dropout mass is deposited non- 
homologously, the dynamic mass to light ratio 
varies with galactic radius, i.e. M(r e /3)/Ls ^ 
Af(480 pc)/ Lb- Even when the dropout pa- 
rameters are scaled homologously from 4472 
[(q , r do , m) = (4, 2 kpc, 1)] to 4472/4 [(q 
= (4, 0.667 kpc, 1)], the resulting M/Lb are not 
equal, i.e. dropout does not scale homologously. 
Some models in Table 1 were performed with no 
additional secondary infall gas following the time 
t* of galaxy formation (NO in column 5). This 



procedure is justified by the strong correlation of 
normalized X-ray sizes r ex /r e with L x /Lb noted 
by Mathews & Brighenti (1998). We explained 
this correlation with models of tidally truncated 
galactic dark halos and (mostly secondary in- 
fallen) halo gas in smaller galaxies and the transfer 
of this material to halos of more massive, cen- 
trally located elliptical galaxies, a process which 
is dynamically expected in galaxy groups where 
elliptical galaxies are formed. The removal of sec- 
ondary infallen halo gas is therefore more relevant 
to 4472/4 and 4472/16 than for 4472 (which may 
have acquired additional halo gas from neighbor- 
ing galaxies); note that this late arriving cosmic 
gas also has a larger influence on M/Lb for the 
smaller galaxies. 

The non-homologous dropout behavior listed in 
Table 1 is illustrated in Figure 2 in which the 4472 

4472/4 -> 4472/16 variation of M/L B (r e /3) 
and M/L B (480 pc) are plotted for each model. 
The local fundamental plane, M/L B oc L B , 
shown with dashed lines, has a slope (3 = 0.24 
that is intermediate between values for the B-band 
FP taken from Jorgensen, Franx, & Kjaergaard 
(1996) and Scodeggio et al. (1998). The dashed 
lines are normalized to M/Lb = 9.2 (from van der 
Marel 1991; long dashed line) and M/L B = 10.3 
(from Faber et al. 1997; short dashed line) re- 
spectively. The asterisks show the constant stellar 
M/ L B = 6 assumed for luminous stars in all three 
galaxy models. Clearly, none of the model pairs 
exhibit a variation of M/Ls(r e /3) as steep as the 
FP relation, but the dropout contribution to the 
FP slope is nevertheless consequential at this ra- 
dius. However, the variation of M/Lb at r = 480 
pc for the q Q = 1 models follows the fundamen- 
tal plane quite accurately over the range of Lb 
between 4472 and 4472/4. Note that the pair of 
variable dropout models with homologously scaled 
parameters [(q , rdo, to) = (4,2 kpc, 1) for 4472 
and (q , r do , m) = (4,0.667 kpc, 1) for 4472/4] 
produce M/Lb (480) that vary in a sense oppo- 
site to the FP trend, but only when all cosmic 
gas from secondary infall is included. Our mass 
to light values for 4472 may be lower limits since 
we have not considered models in which unusually 
massive galaxies receive more than their cosmolog- 
ical allotment of halo material by tidal acquisitions 
from less massive, unmerged group galaxies. 

The constant q = q a = 1 models can explain 
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much of the FP tilt in 4472 -> 4472/4 with or 
without secondary infall. However, when extended 
to 4472/16, values of M/L B (480) clearly exceed 
the fundamental plane for all models. While the 
dropout contribution to M/Ls(480) in the q = 1 
model (without cosmic inflow) decreases with Lb 
over the range of 16 shown in Figure 2 (i.e. vertical 
distance between open squares and asterisks), the 
dropout mass does not disappear entirely for the 
4472/16 model. In fact the assumed M/L B = 6 
for the luminous stars alone (asterisks in Fig. 2) 
also lies above the FP for 4472/16. Within the 
framework of our optically dark dropout assump- 
tion, this indicates that an independent variation 
of M/Lb for the old stars would be required in 
addition to that provided by the dropout. Results 
shown in the upper panel of Figure 2 suggest that 
the FP could be fit over the full 4472 -> 4472/16 
range if we had chosen a somewhat lower value of 
M/L B for the luminous stars, but the X-ray data 
for NGC 4472 become difficult to fit if M/L B < 4. 

However, dynamically significant galactic rota- 
tion is expected for elliptical galaxies having B- 
band luminosities less than that of 4472/4, 2 x 
10 10 Lb,q, which lies near the transition between 
massive, boxy elliptical galaxies with cores and 
low mass, rotationally supported, disky ellipti- 
cal galaxies with power law profiles (Faber et al. 
1997). Rotation can reduce the mass dropout in 
smaller galaxies since the cooled gas will be de- 
posited in a relatively larger galactic volume, re- 
ducing M/L values for 4472/16 in Figure 2. Al- 
though the excess interstellar density in our mod- 
els within 1-3 kpc (Figure 1) can be understood 
in terms of ~ 100 fiG magnetic fields, it could also 
indicate that we have somewhat overestimated the 
gas density in this region and also the mass of 
cooled gas. 

Perhaps the best resolution of the large M/L 
in our models of 4472/16 would be to allow some 
of the cooled gas to form a luminous stellar pop- 
ulation. Using this assumption, we have shown 
that the anomalously low ages for bright elliptical 
galaxies found from observations of the H/3 pho- 
tometric index (e.g., Trager 1997) can be quite 
naturally produced by a young population with 
a Salpeter IMF extending up to ~ 2 M Q but no 
further (Mathews & Brighcnti 1999a). In mod- 
els having luminous dropout stars, about 10 per- 
cent of the optical light is produced by the young 



stars. The stellar disk component characteristic 
of rotating, low luminosity elliptical galaxies (Rix, 
Carollo, & Freeman 1999), can be formed from 
younger cooling dropout stars (Brighcnti & Math- 
ews 1997b), and there is some evidence that these 
galaxies also have enhanced H/3 indices (de Jong 
& Davies 1997). 

In our models we determine the detailed radial 
variation of the mass to light ratio within r e ; this 
variation depends most sensitively on our choice 
of q(r). In general, optical observations of stel- 
lar velocity dispersion profiles cannot determine 
the detailed variation of M/L with radius within 
r e . Only with considerable effort have observers 
been able to detect the presence of dark matter at 
r>r e (e.g., Carollo et al. 1995; Rix et al. 1997), 
although this dark matter is immediately apparent 
from X-ray observations (e.g. Brighenti & Math- 
ews 1997a). The insensitivity of optical observa- 
tions to the detailed radial variation of M/L is 
probably due to the generic stellar orbital bias to- 
ward radial anisotropy in r ^ r e ; this results in 
M/L determinations that more nearly represent 
orbital averages rather than local values. 

For completeness, in the final column of Table 
1 we list the ROSAT X-ray luminosity of each of 
our models. In general our results arc slightly flat- 
ter than the observed trend for luminous elliptical 
galaxies, L x oc L P B with p w 2, in which the X-ray 
emission is interstellar (e.g. Eskridge, Fabbiano 
& Kim 1995). However, we have shown elsewhere 
(Mathews & Brighenti 1998) that L x /Lb corre- 
lates with the size of the X-ray image and that this 
correlation can be explained by the tidal transfer 
of halo material in galaxy groups from subordi- 
nate elliptical galaxies to the dominant central el- 
liptical in the group. Galactic rotation can also 
influence L x /Lb- Because of these environmental 
influences, the values of L x in Table 1 should not 
be regarded as a strong test of a particular model. 

6. FINAL REMARKS AND CONCLU- 
SIONS 

The full extent of the fundamental plane for el- 
liptical galaxies spans a range of Lb about ~ 3 
times greater than the factor of 16 shown in Fig- 
ure 2 (Burstein et al. 1997). Maintaining the FP- 
HP tilt over this larger range would be difficult 
if cooling dropout is an important contributor of 
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dark mass. Nevertheless, the presence of appre- 
ciable dropout mass seems unavoidable. 

The total mass of cooled interstellar gas present 
in elliptical galaxies could be reduced by supernova- 
driven galactic winds or by the formation of stars 
that have since evolved. For our models we assume 
a Type la rate SNu(t) = SNu(t„)(i„/t) p super- 
novae per 100 years per 10 10 Lb with p — 1 and 
SNu(i„) = 0.03 (Brighenti & Mathews 2000a). 
Winds produced by more numerous Type la su- 
pernovae at early times (i.e. p > 1.3) cannot 
be invoked to expel some of the dropout mate- 
rial from galaxies like 4472/4 or 4472/16 since 
these same higher Type la rates in more mas- 
sive galaxies would greatly overproduce iron in 
the current interstellar gas, in violation of X-ray 
observations. Our value SNu(t„) = 0.03 was cho- 
sen to be consistent with the low interstellar iron 
abundances for elliptical galaxies reported in most 
of the recent literature, zp e ^ ism < 0.4 solar (e.g., 
Tanaka et al. 1994). If the abundances are more 
nearly solar (Buote 1999), as is likely, then a pro- 
portionally larger SNu(i„) is indicated; this can 
result in galactic winds in the 4472/16 model at 
late times t ~ t n . Therefore, the total mass of 
interstellar gas that cools can be reduced by SNII- 
driven galactic outflow at early times (when most 
of the dropout occurs) and (less efficiently) by 
SNIa-driven outflow at late times. However, the 
large L x observed in virtually all bright elliptical 
galaxies clearly rules out ongoing galactic winds 
in these galaxies. The current amount of dropout 
mass could also be reduced if the IMF of the 
forming stars was more luminous in the past (i.e. 
more massive stars). While none of these means 
of reducing the dropout mass is securely under- 
stood, some significant mass of cooled gas must 
be present. The large mass of cooled interstellar 
gas expected within r e together with the alarming 
number of possible ways of enhancing or moderat- 
ing the amount of cooled gas makes the coherency 
and tightness of the FP quite baffling. 

Our main objective in presenting these results 
is to emphasize the importance of cooling dropout 
mass to the FP-HP controversy. The results that 
we have shown here are based on the commonly 
held assumption that only optically dark stars of 
very low mass can form in cooling flows. How- 
ever, optically dark dropout by itself appears to 
be an unpromising explanation for the full FP tilt, 



although with some well chosen parameters that 
might be possible. What we see in Figure 2 is an 
unwelcome and quite unaesthetic disturbance of 
the M/L relation, one that involves the uncertain 
physics of star formation and the origin of entropy 
fluctuations in the hot gas. 

Nevertheless, the mass contribution from cooled 
interstellar gas cannot be dismissed by those con- 
cerned with the deviation of the fundamental 
plane from a truly homologous plane. Like other 
effects proposed to explain the FP-HP discrep- 
ancy, if the cooling dropout mass is a dominant 
contributor, the mass dropout must vary in a 
smooth, regular fashion across the FP, propor- 
tional to some power of Lb- Such regularity could 
conceivably result from the similar profiles of in- 
terstellar gas density and temperature observed in 
the central regions of elliptical galaxies (r ^ r e ), 
provided the entropy fluctuations are also simi- 
lar. Our apparent success in explaining the the 
anomalously young ages of elliptical galaxies indi- 
cated by the H/3 index from an optically luminous 
dropout stellar population (Mathews & Brighenti 
1999b) strongly suggests that a better match to 
the FP-HP tilt could be achieved if the optically 
dark assumption used here is relaxed. In this case 
the IMF of the continuously forming stars is likely 
to extend to higher maximum masses in elliptical 
galaxies of lower L b because of the lower interstel- 
lar pressure in low luminosity elliptical galaxies. 

Studies of the evolution of hot gas in elliptical 
galaxies at UC Santa Cruz is supported by grants 
from NASA and the NSF for which we are very 
grateful. In addition FB is supported in part by 
Grant MURST-Cofin 98. 
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Table 1 

COOLING DROPOUT PARAMETERS AND MASS TO LIGHT RATIOS 



GALAXY 


Lb 
(10 10 i B ,o) 




r d o a 
(kpc) 


77J a 


COSMIC b 
GAS? 


M (r e \ 
L B <■ 3 / 


£(480 pc) 


L x c 

(10 40 ergs s- 1 ) 


4472/16: 






















0.493 









NO 


9.32 




8.57 


0.10 




0.493 


1 


oo 




NO 


7.29 




7.13 


0.10 




493 


4 


317 

U.ul 1 




NO 


8 25 




8.17 


25 


4472/4: 






















1.97 









YES 


11.76 




19.27 


9.35 




1.97 









NO 


8.97 




12.37 


1.56 




1.97 


1 


OO 




YES 


7.95 




8.41 


8.00 




1.97 


1 


oo 




NO 


7.40 




7.65 


1.52 




1.97 


4 


0.667 


1 


YES 


11.05 




10.91 


14.57 




1.97 


4 


0.667 


1 


NO 


8.65 




8.54 


2.57 




1.97 


4 


2 


1 


YES 


8.16 




7.27 


10.78 




1.97 


4 


oo 




YES 


7.03 




6.79 


6.00 


4472: 






















7.89 









YES 


10.47 




32.29 


20.08 




7.89 









NO 


9.97 




28.53 


8.59 




7.89 


1 


oo 




YES 


8.83 




11.40 


31.56 




7.89 


1 


oo 




NO 


8.64 




10.96 


12.56 




7.89 


4 


2 


1 


YES 


9.74 




7.99 


28.11 




7.89 


4 


2 


1 


NO 


9.37 




7.88 


13.73 




7.89 


4 


0.8 


2 


YES 


10.47 




13.08 


36.44 




7.89 


4 


oo 




YES 


7.88 




7.20 


39.82 



a Whcn q(r) = q a is constant, the exponential factor in equation (3) is set to unity, as if rdo = oo and m is 
undefined. 

b YES (NO) indicates that continued secondary flow of cosmic gas into the galactic halo is (not) included. 

c ROSAT band X-ray luminosity (including dropout emission) within 150 kpc for NGC 4472, 50 kpc for 
4472/4 and 16.6 kpc for 4472/16 
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Fig. 1. — Top panel: Comparison of radial in- 
terstellar density variation for two hydrodynamic 
models at time t n = 13 Gyrs with observed den- 
sity in NGC 4472. Filled circles are Einstein data 
(Trinchieri, Fabbiano, & Canizares 1986), open 
circles are from ROSAT (Irwin & Sarazin 1996). 
The light solid and dashed lines refer to the back- 
ground (uncooled) gas density in the q a = 1 and 
(<Zo, rdo, to) = (4, 2 kpc, 1) models respectively. 
The heavy lines are the apparent space variation 
of the gas density including emission from gas that 
is cooling locally. Bottom panel: Comparison of 
the interstellar temperature variation with radius 
computed in the same two models with ROSAT 
temperature data (Irwin & Sarazin 1996). The 
light solid and dashed lines refer to the computed 
temperature variation in the uncooled background 
flow with physical radius. The heavy solid and 
dashed lines show the apparent variation of tem- 
perature projected on the sky and corrected for 
emission from mass dropout; these curves should 
be compared with the observations. 

Fig. 2. — Variation of the mass to light ratio 
M/Lb with luminosity Lb (both in solar units) 
at t n = 13 Gyrs. M/Lb includes both luminous 
stars and interstellar dropout for NGC 4472 (on 
the right) and two homologously scaled galaxies 
having 1/4 and 1/16 the mass of NGC 4472 (cen- 
ter and left respectively). Top panel: M/Lb eval- 
uated within fixed radius r = 480 pc for all galax- 
ies; Bottom panel: M/Lb within r = r e /3 for all 
galaxies. Asterisks show M/Lb = 6 assumed for 
the luminous stars. Dashed lines show the local 
the fundamental plane, M/Lb °c Lb ' 24 , normal- 
ized at NGC 4472 with M/L B = 9.2 (long dashed 
line) or M/Lb = 10.3 (short dashed line). Circles 
refer to no-dropout models (q D = 0), squares to 
uniform q Q = 1 models, triangles to uniform q D = 4 
models and hexagons to variable dropout models 
with (q D , rdo, m) = (4, 2 kpc, 1) for 4472 and 
(q Q , r do , m) = (4,0.667 kpc, 1) for 4472/4 (up- 
per filled hexagon) or (q Ql rd 0l m) = (4,2 kpc, 1) 
for 4472/4 (lower filled hexagon). Gas from sec- 
ondary cosmic infall is included in solutions with 
filled symbols and this gas is excluded in solutions 
with open symbols. Where possible the symbols 
are labeled with (q D , r^o, m, Y/N); Y/N refers 
to the presence or lack of secondary infall. 
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